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ABSTRACT 

We present Goddard High Resolution Spectrograph observations of the M8 Ve star VB 10 ( = GI 752B) 
located very near the end of the stellar main sequence, and its dM3.5 binary companion G1 752A. These 
coeval stars provide a test bed for studying whether the outer atmospheres of stars respond to changes in 
interna 1 structure as stars become fully convective near mass 0.3 M Q (about spectral type M5), where the 
nature of the stellar magnetic dynamo presumably changes, and near the transition from red to brown dwarfs 
ear mass 0.08 (about spectral type M9), when hydrogen burning ceases at the end of the main sequence 
We obtain upper limits for the quiescent emission of VB 10 but observe a transition region spectrum during a 
large flare, which indicates that some type of magnetic dynamo must be present. Two indirect lines of 
evidence scaling from the observed X-ray emission and scaling from a time-resolved flare on AD Leo— 
suggest that the fraction of the stellar bolometric luminosity that heats the transition region of VB 10 outside 
of obvious flares is comparable to, or larger than, that for G1 752A. This suggests an increase in the magnetic 
heating rates, as measured by L, lne /L bol ratios, across the radiative/convective core boundary and as stars 
approach the red/brown dwarf boundary. These results provide new constraints for dynamo models and 
models of coronal and transition-region heating in late-type stars. 

Subject headings: binaries: visual — stars: activity — stars: flare — stars: individual (VB 10) — 
stars: interiors — stars: late-type — X-rays: stars 


1. INTRODUCTION 

Two fundamental changes in the internal structure of stars 
occur near the end of the hydrogen-burning main sequence. 
With decreasing stellar mass and T eff , the radiative core 
shrinks until, for masses below about 0.3 Af 0 (corresponding 
to = 3400 K or spectral type M5), stars are fully convective 
(see Dorman, Nelson, & Chau 1989). This structural change is 
important because the shell-type a-co magnetic dynamo 
described by Parker (1979) and others, which is thought to 
amplify magnetic fields in solar-type stars, requires radial dif- 
ferential rotation. The important result from helioseismology 
that the solar convective zone has nearly constant angular 
velocity with depth led Parker (1993) and others to place the 
radial differential rotation that creates toroidal magnetic fields 
(the “cu-effect ”) not in the convective zone but rather in a layer 
of about 0.03 R Q thickness at the top of the radiative zone. It is 
not presently known whether the a-co dynamo or another type 
of dynamo — for example, the “distributive” dynamo discussed 
in general terms by Rosner (1980)— can operate in fully con- 
vective stars. However, Durney, De Young, & Roxburgh (1993) 
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showed that fully convective stars can generate turbulent but 
not global magnetic fields. 

As one proceeds further down the main sequence to 0.076 < 
< 0.080 (near spectral type M9), depending on the 
assumed equation of state, solar composition stars can no 
longer support stable hydrogen burning (Dorman et al. 1989; 
Burrows, Hubbard, & Lunine 1994). Such stars are commonly 
called brown dwarfs, although a candidate brown dwarf has 
yet to be identified with certainty. 

Although it is not presently possible to study these changes 
in internal structure directly, the properties of stellar magnetic 
fields and the magnetic heating of the outer atmospheric layers 
may be proxy indicators of structural changes occurring deep 
inside the stars. To test this hypothesis, Barbera et al. (1993) 
analyzed X-ray observations of all K and M stars within 25 pc 
that were observed by the Einstein Observatory imaging pro- 
portional counter (IPC). Although this data set is biased by the 
requirements of the individual pointings and there are only a 
few late M dwarf detections and upper limits, Barbera et al. 
found a definite trend of lower X-ray luminosities for the lower 
mass fully convective stars compared to the higher mass M 
dwarfs with radiative cores. This result can arise simply from 
the rapid decrease in stellar radius with decreasing mass com- 
bined with a roughly constant X-ray surface flux (see, e.g., 
Fleming et al. 1993; Fleming, Schmitt, & Giampapa 1995)! 
Barbera et al. also claimed to find a drop in the mean X-ray 
surface flux across the boundary between the early M dwarfs 
with radiative cores and the fully convective lower mass stars, 
but the evidence that they presented is only marginal. 

Fleming et al. (1995) and Schmitt, Fleming, & Giampapa 
(1995) came to a different conclusion in their analysis of the 
ROSAT all-sky survey and pointed observations of all known 
K and M stars within 7 pc of the Sun. Using L^/L^i as a 
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measure of coronal heating efficiency, they found no difference 
between typical values for the early and late M-type stars, 
which indicates that coronal magnetic heating rates are insen- 
sitive to changes in the internal structures of the stars, at least 
at the level detectable by ROSAT. Moreover, the increased 
sensitivity of ROSAT yielded more detections among lower 
mass stars in contrast to the paucity of detections available in 
the Einstein data set. ROSAT did not detect any stars later 
than M7, however, so they could not extend their results to the 
lowest mass stars. Fleming et al. (1993) came to similar conclu- 
sions in their analysis of all known dwarfs later than M5 in the 
ROSAT survey and pointed observations. 

2. THE TARGET STARS VB 10 G1 752A 

We have obtained ultraviolet spectra of the dM3.5 and M8 
Ve stars in the G1 752 binary system in order to study possible 
differences in the amount of plasma heated to 10 5 K in the 
outer atmospheres of two stars with the same age but different 
internal structure. G1 752A has a radiative core, while G1 752B 
(also known as VB 10) is fully convective and lies very close to 
the end of the hydrogen-burning main sequence. G1 752B is the 
faintest (m v = 17.38) and coolest dwarf star yet observed with 
the Goddard High Resolution Spectrograph (GHRS). 

Liebert & Probst (1987) have reviewed the properties of low- 
mass M dwarfs. Using optical and infrared photometry, Berri- 
man, Reid, & Leggett (1992) derived an effective temperature 
and luminosity of T cff = 3240 K and log (LJL 0 ) = — 1.49 for 
G1 752A. Kirkpatrick et al. (1993) derived a temperature and 
luminosity of T eff = 2875 K and log {LJL 0 ) = —3.31 for VB 
10, by comparison of the observed spectrum with model atmo- 
spheres. Tinney, Mould, & Reid (1993) used infrared photo- 
metry to determine the absolute bolometric magnitude of VB 
10, M bol = 13.13, corresponding to log (LJL Q ) = -3.36. 
Newer model atmosphere calculations with a better represen- 
tation of the H 2 0 opacity led Allard et al. (1994) to revise the 
temperature for VB 10 down to T €ff — 2800 K. Most recently, 
Allard & Hauschildt (1995) obtained an even lower tem- 
perature, T cff = 2600 K, by matching the observed near-IR 
spectrum with model atmospheres using new TiO opacities. 

We now derive stellar radii from the relation 

5 log (RJRo) = M bolt0 - - 10 log ( W^.o) • U) 

Adopting M bo i i0 = 4.72 (based on Af K0 = 4.83 [Bessell 1991] 
and BC g = — 0.1 1, which is uncertain), Af bol V B to ~ 13.13, and 
il/boi.Gi 752 A = 8.54 (Kirkpatrick et al. 1993), we derive radii for 
G1 752A and VB 10 of 0.546 and 0.102 R 0 , respectively. The 
sensitivity of the derived stellar radii to the assumed value of 
T« f is indicated by R V bio decreasing from 0.102 to 0.080 R 0 
when T c(f is raised from 2600 to 2800 K. 

The mass of VB 10 is not known directly because it is a 
member of a long-period binary system. Henry & McCarthy 
(1993) list it among low-luminosity dwarfs with colors that 
suggest, following their empirical mass-luminosity relations, a 
mass in the range of 0.08-0.09 Af 0 . The theoretical mass- 
luminosity relations for different assumed equations of state 
(see Dorman et al. 1989) also indicate a VB 10 mass of about 
0.09 Af 0 . 

We derive the ratio of surface flux to observed flux, F// = 
2.17 x 10 17 for G1 752A and 6.23 x 10 18 for VB 10, from the 
stellar radii and the 5.74 + 0.16 pc distance to G1 752 
(Harrington, Kallarakal, & Dahn 1983). The stellar properties 
of Gl 752A and VB 10 are summarized in Table 1. For com- 


TABLE l 


Summary of Adopted Stellar Parameters 


Parameter 

AU Mic 

Gl 752A 
(Wolf 1055) 

Gl 752B 
(VB 10) 

R.A.(2000) 


\9 b \6 m 55'9 

19 h 16 m 58?25 

Decl.(2000) 


4- 5°10'30'.'9 

+ 5 C1 09'24'.'3 

Spectral type 

dMOe 

dM3.5 

M8 Ve 

T cff (K) 

3770 

3240 

2600 



7.89 

8.54 

13.13 

My 


10.31 

18.57 

log(V^o) 

-1.23 

-1.53 

-3.36 

RJRq 

0.544 

0.546 

0.102 

Distance (pc) 

9.35 

5.74 

5.74 

F//(10 17 ) 

5.81 

2.17 

62.3 

M JM q 

0.54 

0.46 

*0.09 

log L x (ergs s ‘) 

29.8 

27.25 

26.34 

1°8 (^x/^bol) 

-2.47 

-4.80 

-3.88 


parison we also list the properties of AU Mic (Linsky & Wood 
1994), which is among the most active of the M dwarfs with 
very bright transition-region ultraviolet emission lines (Linsky 
et al, 1982) and the largest X-ray luminosity of any M dwarf 
within 10 pc. 

The inferred radius for VB 10 is very similar to that of 
Jupiter (R Jup = 0.1024 R 0 ; Allen 1973, p. 140). This is not 
unreasonable, however, because the theoretical calculations of 
Burrows, Hubbard, & Lunine (1989) predict that stars near the 
lower hydrogen-burning mass limit (M » 0.08 Af 0 ) should 
have radii of about 0.08 R 0 (see also Fig. 1 in Burrows & 
Liebert 1993). VB 10 is certainly near this limit. Kirkpatrick & 
McCarthy (1994) raise the possibility that stars with spectral 
types M7 and later are, in fact, below this limit; thus, VB 10 
could conceivably be a brown dwarf. 

X-ray emission from Gl 752A was detected by the Einstein 
IPC at a level of log L x = 27.25 (Barbera et al. 1993) and by the 
ROSAT Position Sensitive Proportional Counter (PSPC) at a 
level of log L x — 26.71 (Schmitt et al. 1995). Since (for the IPC 
data) log (Ex/Lfoi) = -4.80, nearly 2 orders of magnitude 
below the “saturated” heating level [log (LJL hol ) ~ 10“ 3 ; see, 
e.g., Fleming et al. 1995, their Fig. la , and Gagne, Caillault, & 
Stauffer 1995], Gl 752A is relatively inactive. Fleming et al. 
(1993) provided an upper limit of log L x < 25.54 (L x /L bol < 
— 4.68) for VB 10 based on a long pointed ROSAT PSPC 
observation, but Barbera et al. (1993) cited an Einstein HRI 
detection corresponding to log L x = 26.34 in the IPC band. 
The corresponding value of log (L x /L bol ) = —3.88 is an order 
of magnitude larger than that of Gl 752A, which suggests an 
increase rather than a decrease in the coronal heating efficiency 
across the radiative/nonradiative boundary. The X-ray detec- 
tion of VB 10 with the HRI, but not the IPC or the PSPC, 
suggests variability or flaring. T. A. Fleming (1994, private 
communication) has looked at the underexposed 9.1 ks HRI 
observation obtained in two segments (due to observing 
constraints) over an 8 h period, but found no evidence for flares 
on shorter timescales in the arrival time of the seven detected 
photons. Clearly additional observations are needed to 
confirm the X-ray luminosity of VB 10 and to study its likely 
X-ray flares. 

The amount of Ha absorption or emission indicates the 
chromospheric heating rate, which is correlated with stellar 
kinematic properties and thus age (e.g., Giampapa & Liebert 
1986). Gl 752A shows a modest Ha absorption equivalent 
width of 0.49 A (Panagi, Byrne, & Houdebine 1991), which 
indicates a relatively low amount of heating for stars of its 
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has a very different character with prominent emission lines of 
C ii. Si iv, and C iv. This is clearly a flare spectrum. Although 
VB 10 is known to flare in the optical (see, e.g., Herbig 1956) 
and X-ray flares have been detected in stars with masses nearly 
as small as VB 10 (e.g., the flare on VB 8 observed by Taglia- 
ferri, Doyle, & Giommi 1990), this is the first detection of 
flaring high-temperature gas in the ultraviolet from such a 
low-mass star. Upper limits for the line fluxes and surface 
fluxes in the quiescent spectrum and the detected fluxes in the 
flare spectrum are listed in Table 3. 

Because our G1 752 spectra are quarter-stepped, each 
readout actually consists of four separate sequential exposures, 
each of 81 s duration. We can therefore study the time evolu- 
tion of the VB 10 flare during readout 1 1 by measuring from 
the raw data the number of C iv counts detected during each of 
the four 81 s exposures of the readout. The number of C iv 
counts observed during each sequential exposure are 33 + 6, 
38 ± 6, 31 ±6, and 24 ± 5. These numbers suggest that the 
C iv flux was decreasing by the end of readout 11. Because 
there is no trace of the flare present in readout 10, the flare 
must have started during the 31 s dead time between readouts 
10 and 1 1, or during the first exposure of readout 11. The time 
evolution of the flare was not fully sampled, but the decrease in 
C iv flux during the last two 81 s elements of the exposure 
suggests a 1/e folding time of At « 5 minutes. The rapid rise 
and decay of the flare is similar to those of solar impulsive-type 
flares. 

For post-COSTAR G140L spectra, the expected spectral 
resolution is about 0.65 A (Soderblom 1994). The widths of the 
emission lines in the G1 752 A and VB 10 spectra are consistent 
with this resolution, with the exception of the C iv lines 
observed during the flare on VB 10. The widths of these lines, 
as measured by Gaussian fits, are about 1.3 A, indicating that 
the intrinsic line widths during the flare are more than 1.0 A, 
corresponding to non-thermal velocities of roughly 80 km s~ 1 
or more. These velocities are comparable to velocities inferred 
from solar transition-region lines during solar impulsive flares 
(Cheng 1990). We did not observe any large redshifts, which 
presumably indicate rapidly downflowing matter, as Book- 
binder, Walter, & Brown (1992) found in their GHRS spectra 
of the 1991 May 9 flare on AD Leo. 

We estimate the emission measure of the flaring 10 5 K 
plasma from the C iv surface flux using the relation EM = 
Fdv/2.056 x 10“ 23 = 1.1 x 10 28 cm 5 (see Brown et al. 
1984). The corresponding volume emission measure is 

VEM = 4/rR 2 EM = 6.8 x 10 26 (^/10‘ l )~ 2 cm ” 3 . (2) 

The corresponding linear dimension is 

x = (VEM)‘ 1/3 = 8.8 x 10 3 (n e /10 u r 2/3 km 

= 0.13R J> (n e /10 11 )' 2/3 . (3) 

Thus, the 10 5 K flaring plasma covered a large fraction of the 
star if n e % 10 11 cm” 3 , but only a small fraction if n e > 10 12 
cm” 3 . We note that the distance traveled by a sound wave in a 
10 5 K plasma in 5 minutes, the estimated \/e decay time of the 
flare, is about 0.1R*. 

4.2. T he Nonflare Spectrum of VB JO 

Table 3 contains the line surface fluxes (flux per unit area of 
the stellar photosphere) obtained using the F/f ratios in Table 
1. In Table 4 we compare the ratios of the line luminosity to the 
bolometric luminosity, L^JL^, for the G1 752 stars to those 


for the quiet Sun, AU Mic, and the fully convective dM5.5e 
star, Proxima Centauri, observed outside of obvious flares 
using the IUE fluxes tabulated by Linsky et al. (1982). We 
intercompare these stars using because this ratio 

measures the heating efficiency where the line is formed, but, 
unlike the surface flux, this ratio is independent of distance and 
stellar radius, which is decreasing rapidly for the very late M 
stars. 

Despite the 54 minute integration for the nonflare VB 10 
spectrum, the 2 o upper limits on the L^JL^i ratios are larger, 
typically by an order of magnitude, than the corresponding 
detected values for the companion star G1 752A. The large 
upper limits are due in part to the factor of 30 decrease in the 
square of the radius of VB 10 compared to that of G1 752A. If 
the value of Lcy/L^ for VB 10 were the same as for G1 752A, 
for example, then we would have to integrate for 40 hr to detect 
the C iv line in VB 10 considering only Poisson statistics in the 
line flux; but, in reality, we would have to integrate for a much 
longer time, since the noise in the background dominates the 
error budget. In order to speculate whether the nonflare 
plasma heating rates at transition-region temperatures change 
across the radiative/convective core boundary or as stars 
approach the red/brown dwarf boundary, we consider two 
indirect methods for inferring the nonflare L line /L bol ratios. 

First, we note that the Einstein HRI value of LyJL M for 
VB 10 listed in Table 4 is a factor of ~ 1.8 smaller than for the 
quiescent Prox Cen, and the upper limits in the nonflare VB 10 
^line/^boi ratios are about a factor of 5 larger than the corre- 
sponding measured values for Prox Cen. Because the X-ray 
and transition-region line emission are observed to be related 
by a power law of the form L x ^ (Ayres, Marstad, & 
Linsky 1981), we estimate that the nonflare values for the VB 
10 transition-region lines should lie a factor of 1.5 below the 
values for Prox Cen and thus a factor of 5 x 1.5 = 7.5 below 
the VB 10 upper limits in Table 4. This line of reasoning leads 
to the quiescent values of L, ine /L bol being a factor of 10/7.5 = 
1.3 larger in VB 10 than in G! 752A. Since the value of L x for 
VB 10 is uncertain and may include flares, this argument by 
itself does not carry much weight. 

The flare on VB 10 provides data for a second way of esti- 
mating L linc /L bo , for VB 10 outside of flares. In Table 5, we 
compare flare peak C iv luminosities and total flare lumi- 
nosities L nare emitted by plasma hotter than 10 4 K for three 
late M stars, using the empirical relation between the C iv 
emission and the total radiative losses in the transition region 
and corona proposed by Doyle (1989) for M dwarfs. The C iv 
fluxes for the AD Leo flare observed with the GHRS on 1991 
May 9 (Bookbinder et al. 1992) refer to the flare peak emission, 
as is the case for the VB 10 flare, because the GHRS provides 
sufficient time resolution to measure the peak emission. On the 
other hand, the C iv fluxes measured from IUE spectra for the 
Prox Cen flare (Haisch et al. 1983) and the giant 1985 April 12 
flare on AD Leo (Hawley & Pettersen 1991) refer to 60 minute 
and 15 minute time averages, respectively, which could be sig- 
nificantly smaller than the true flare peak fluxes. For this 
reason, we consider the IUE fluxes and the derived lumi- 
nosities for these two flares to be lower limits to the flare peak 
values. 

We take as our guide the 1991 May 9 flare on AD Leo, 
because, like the VB 10 flare, the peak value of f c , v is measured 
and the flare is of short duration. Although the ratio of the 
peak C iv flux to the quiescent value, /£j V //civ> must surely 
vary from flare to flare, the ratio of 90 observed in the flare on 
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TABLE 4 

Values of log { L Vl JL ^ 

A mult G1 752A VB 10 VB 10 Quiet Sun AU Mic Proxima Centauri 


Ion (A) (dM3.5) (M8 Ve) (flare) (G3 V) (dMOe) (dM5.5e) 


X-rays -4.80 -3.88 ... -6.30 -2.47 -3.62 

Nv 1240 -7.10 ... ... -8.60 -6.03 -6.58 

Civ 1549 -6.51 <-5.54 -4.08 -7.59 -5.46 -6.25 

Si iv 1398 -7.25 <-5.84 -4.89 -7.85 -5.88 

Cn 1335 -6.61 <-5.84 -4.97 -7.58 -5.90 -6.53 

Hen 1640 -6.78 ... ... -8.33 -5.78 -6.75 

Oi 1304 ... <-5.67 <-5.28 -7.66 -5.96 -6.75 


AD Leo may be representative of large flares. The data in 
Table 3 show that the ratio must exceed 40 for the VB 10 flare, 
since the observed C iv 1548 A flux is 40 time larger during the 
flare than the nonflare upper limit. Assuming the value of 90 
for the ratio, we estimate that the quiescent value for VB 10 is 
fcw = 4.0 x 1(T 16 ergs cm“ 2 s _1 , and thus the quiescent 
Ldv/Lboi ratio is about two-thirds of the measured upper limit. 
This line of reasoning suggests that the quiescent L clv /L boj 
ratio is a factor of 7 larger for VB 10 than for G1 752. 

Although falling short of proof, these two indirect but inde- 
pendent lines of evidence — scaling from the observed X-ray 
emission and scaling from a time-resolved flare on AD Leo — 
both indicate that the fraction of the stellar bolometric lumi- 
nosity that heats the transition region of VB 10 outside of 
obvious flares is comparable to or larger than that in G1 752A. 
Because the two stars have the same age, we cannot explain 
this result as an age effect. Instead, our analysis suggests an 
increase in the magnetic heating rates across the radiative/ 
convective core boundary and as stars approach the red/brown 
dwarf boundary, although the degree of this increase is uncer- 
tain. Another way of looking at this is that in order for the 
magnetic heating rates across the radiative/convective core 
boundary to decrease, the /civ//civ ratio must exceed 630 and 
the X-ray scaling law of Ayres et al. (1981) must be invalid for 
late M dwarfs. There is no evidence to support either of these 
arguments. 

5. CONCLUSIONS 

We have observed ultraviolet emission lines indicating 
plasma heated to 10 5 K during a flare on VB 10, the lowest 
mass star showing evidence for transient heating to these tem- 
peratures. The detection of this flare, together with the pre- 
vious optical observation of transient heating of 10 4 K 
chromospheric material, shows that even the lowest mass stars 
at the end of the hydrogen-burning main sequence have strong 
but unstable magnetic fields, and that they must have some 
form of dynamo to replenish the annihilated surface magnetic 
fields. Although we obtain only upper limits for the VB 10 


emission-line fluxes outside of the flare, two indirect lines of 
evidence allow us to estimate crudely what the line fluxes may 
be. These estimates suggest either constant or increasing mag- 
netic heating rates as measured by Ly^JL ^ , ratios across the 
radiative/convective core boundary and as stars approach the 
red/brown dwarf boundary. The X-ray detection of VB 10 by 
the Einstein HRI is consistent with this conclusion, although 
flares could contribute significantly to this detection, but the 
ROSAT upper limit for VB 10 is only at the same L^JL^ level 
as Gl 752A. 

The nature of the nonradiative heating of the outer atmo- 
sphere of VB 10 appears to be characterized by a substantial 
transient component, which we attribute to flare activity. The 
dynamo theory and theories for the electrodynamic heating of 
stellar atmospheres must now account for both the occurence 
of transient flarelike activity with efficient heating rates and the 
inferred increase in the steady quiescent heating rates com- 
pared to those of earlier M dwarfs of the same age and chemi- 
cal composition. As it is difficult to draw definitive conclusions 
from observations of only one system for a limited time, further 
Hubble observations are needed. 

Because the standard a-co type dynamo model apparently 
requires a radiative core to operate, we urge the development 
of detailed models for a magnetic dynamo that can operate in 
fully convective stars. This paper suggests that the dynamo 
may have to be capable of converting a larger fraction of the 
stellar bolometric luminosity into heat for plasma at 10 5 K 
than the standard ot-co type dynamo can in early M stars of the 
same age. 

This work is supported by NASA Interagency Transfer 
S-56500-D to the National Institute of Standards and Tech- 
nology. We thank Richard Robinson for confirming the geo- 
coronal nature of the O i emission observed in the VB 10 
spectrum, Thomas Fleming for bringing the X-ray detection of 
VB 10 to our attention, and France Allard for sending to us her 
new models for VB 10. We also thank James Liebert and Marc 
Gagne for their input. 


TABLE 5 


Flare Peak Luminosities for Late M Dwarfs 


Parameter 

VB 10 

(1994 Oct 12) 

Proxima Centauri 
(1980 Aug 20) 

AD Leo 
(1991 May 9) 

AD Leo 
(1985 Apr 12) 

tog 

30.22 

30.78 

32.00 

32.00 

f F 

J c IV 

3.56 x 10" 14 

>9 x 10 -13 

2.0 x 10 11 

>1.9 x 10" 11 

L C1V 

1.4 x 10 26 

>2 x 10 26 

5.7 x 10 28 

>5 x 10 28 

L c iv/Eboi 

9.0 x 10' 5 

>3 x 10" 5 

5.7 x 10‘ 4 

>5 x 10~ 4 

^flare 

3.9 x 10 31 

>5.5 x 10 31 

2.6 x 10 34 

>2 x 10 34 

Eflare/Eboj 

25 

>9 

260 

>250 

At (minutes) 

*5 

<120 

1 

*60 

/£,v//g,v 

>40 

>7 

90 

>44 
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